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SUMMARY 

Current formulations of the generation of aerodynamic sound 
by turbulence all require statistical Information with regard to 
the turbulent flowfield. Second-order closure turbulent modeling 
in its present form is shown to provide some of this statistical 
Information. 

The Ribner 1 formulation of the generation of aerodynamic 
sound is coupled with predictions of second-order velocity corre- 
lations and integral scale to estimate the sound radiated from 
several complicated jet flows. In particular, it Is shown that 
the sound radiated from a cold swirling jet is greater than from 
its nonswirling equal thrust counterpart. 

The noise radiated from the flowfield of a multitube suppres- 
sor is estimated and compared with an equal thrust equal diameter 
Gaussian jet. It is shown that the multitube concept Is Indeed 
quieter. 


I . INTRODUCTION 


This report describes an attempt to develop the technology 
required to predict the aerodynamic noise radiated from general 
turbulent flowfields. It is but a first attempt at an admittedly 
very difficult problem. The difficulties are associated with our 
quite limited description of turbulence. Nevertheless, It is shown 
that advanced turbulent transport theory, namely second-order 
closure turbulence theory, does Indeed provide new Insight into 
the aerodynamic sound generation problem. It is also suggested 
that second-order closure modeling does provide a framework or 
backbone upon which can be built a more complete turbulence theory, 
which is directly suited to be coupled with accepted formulations 
of the aerodynamic generation problem. The most popular of these 
are the theories of Ribner, 1 Lighthill, 2 Lilley 3 and most recently 
that of Yates. 4 

What is offered by second-order closure turbulence theory is 
a method of computing certain statistical turbulent quantities 


under fairly general circumstances. A deficiency, in our opinion, 
of several models currently being used to estimate the sound 
radiated from a turbulent jet, is that the model utilizes turbulent 
assumptions or methodology which are only applicable to one fluid 
geometry, namely a jet. The introduction of other jets, heat or 
swirl, may result in the modification of the turbulence and, there- 
by, make predictions from these models of limited usefulness. It 
is for this reason that a truly uniform and general approach to 
computation of turbulent transport and resulting generation of 
aerodynamic sound is sought. 

This report is organized as follows. In Section II, we argue 
why at least a second-order closure turbulence theory Is needed 
to compute turbulent quantities to be utilized to estimate the 
aerodynamic generation of sound. In Section III, we show some pre- 
dictions of turbulent properties made by second-order closure theory 
and comparison with data. A simple acoustic model suggested by 
Ribner 1 is reviewed in Section IV, and modifications and assumptions 
required to utilize predictions made from second-order closure 
theory are detailed. In Section V, predictions of the noise radi- 
ated in cold jets are given. Included here are predictions for 
jets which have swirl as well as a multitube jet configuration. 

The effect of heat and swirl on turbulent jet structures is inves- 
tigated in Section VI. Finally in Section VII, conclusions are 
offered. 

NOMENCLATURE 

a Q reference speed of sound 

A axial parameter, see Eq. ( 1 4 ) 

b turbulent model constant 

D jet diameter 

N swirl parameter, see Eq. (14) 

P acoustic power radiated in the <J> direction per unit solid 

angle 


q r.m.s. turbulence level, yu^q^ 

r,0,z circular cylindrical coordinate system 
t time 

velocity components In a Cartesian coordinate system 
u,v,w velocity components in a r,0,z coordinate system 

r p radius of primary jet 


2 



v turbulent model constant 

c 

x.j, Cartesian coordinates 

a swirl angle 

p density 

tt pressure 

v kinematic viscosity 

A turbulent integral scale parameter 


II. SECOND-ORDER CLOSURE TURBULENT TRANSPORT 
THEORY (INCOMPRESSIBLE PLOW) 

Second-order turbulent closure theory is a model by which_one 
can compute the ensemble average mean components of velocity U-^ , 
as well as the second-order turbulent velocity correlations u^u^ . 
The derivation of the second-order closure model is straightforward 
and is detailed briefly below. 


The dynamics of an incompressible fluid are governed by the 
continuity equation and the momentum equations 


3U, 
3x . 


= 0 


( 1 ) 


!Zi + u = 

3t J 3Xj p 3x^ 3 Xj3Xj 


( 2 ) 


Assuming the velocities and pressure can be split as the sum of 
an ensemble average mean and a fluctuation 


u i - u i + u i 


IT = IT + p ' 


( 3 ) 


there results upon substitution into Eqs. (1) and (2) and ensemble 
averaging 
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equations for the mean velocity and pressure. The momentum equa- 
tions, however, contain unkno wn se cond-order velocity correlations 
known as Reynolds stresses -ujM . Solution of Eq s. (4) and (5) 
cannot be attempted until a specification of u£uj is made. The 
essence of first-order clo sure , K theory or eddy viscosity theory 
is an attempt to specify uj^uj as a function of mean flow gradi- 
ents and mean flow length scales. 


Second-order closure theory derives independent equations for 
these second-order velocity correlations. The procedure for doing 
this is as follows. The splitting of variables into mean and 
fluctuation is substituted into Eq. (2), and Eq. (5) is subtracted 
from this result. The remainder is then multiplied by u- and 
ensemble averaged. The above procedure is then repeated inter- 
changing i for j and then it can be shown that the second- 
order velocity correlations must satisfy 
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velocity and pressure 
diffusion 


auM 
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3 u'u' 

v — M- 


3x 
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viscous 

diffusion 


+ !a.' 

P 1 3x 1/ 


pressure strain 
correlation 


- 2v 


3u ! 3u! 


dissipation 


( 6 ) 


Note now that there appears third-order and some second- 
order correlations which themselves need independent equations. In 
principle third-order correlation equations could be developed but 
the spirit of second-order closure is to model these third-order 
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and second-order correlations in terms of known second-order 
correlations. The success of the second-order theory is dependent 
upon two assumptions: 

1. The Reynolds stress equations are not extremely sensi- 
tive to inaccuracies introduced by the replacement of 
exact terms with modeled ones. 

2. Reasonable physically motivated models can be developed 
for the terms which must be replaced. 

These assumption are apparently Justified, since it has been 
demonstrated over the years that the second-order closure model 
has predicted turbulent transport in a variety of high Reynolds 
number flowfields. The details of the modeling as well as compar- 
isons of the turbulent model with data may be found in Ref. 5. 
Here, we simply give the terms to be modeled and the terms which 
replace them. 


Pressure - velocity diffusion 



u k u l u J 



where v has been determined to be 0.3, 
the turbulent macroscale or integral scale 


c 3x. 


qA 
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q = u£ui 
variable . 


i J 


3x k 

) (7) 

and 

A is 


Pressure strain correlation 





( 8 ) 


where this term contains the mechanisms by which the energy com- 
ponents attempt to redistribute energy between themselves. 


Dissipation 


3u’ 3u. ,3 

\) js j — £ _ 0. 

3x k 3x k " °ij 3A (9) 

where b has been determined to be equal to 0.125. The dissipa- 
tion here is for high Reynolds numbers and is assumed to be iso- 
tropic. Upon substitution of these modeled terms the equations 
governing the dynamics of the second-order correlations become 
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( 10 ) 


The turbulent Integral scale parameter A is determined from 

3A , 3A _ , o A — — 8U i . n r „ 

3t U k 3x k ' " 1-8 q 2 u i u k 3x k 0,6 q 


. „ 3 / . 3A 

c 3x k l 3x k 
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The details of the derivation of this relationship is given in 
Ref. 5. 


Equations (4), (5), (10) and (11) along with appropriate 
bound ary c onditions are then solved numerically for the quantities 
u i * u-^Uj and A . 

Before proceeding to obtain solutions for turbulent flowfields 
relevant to the aerodynamic sound generation problem, we can show 
that second-order closure turbulent modeling does contain physics 
that might be lost in K or eddy viscosity theories. We do this 
by asking if second-order closure can be simplified under some 
limiting process to define an eddy viscosity model. The conditions 
under which an eddy viscosity can be expected to be valid have been 
argued by Donaldson. 6 Briefly, eddy viscosity concepts are re- 
stricted to high Reynolds number flows, where second-order correl- 
ations can track their equilibrium values, and the effects of dif- 
fusion of correlations can be neglected. If we now consider a 
nonspreading cylindrical jet with swirl, where the coordina tes a re 
r,0,z_ wi th ve locity components u,v,w , with U = 0 , V = V(r ) 
and W = W(r ) and take the high Reynolds number, nondiffusive 
superequilibrium limit of Eq. (10), there results 
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4uvV - jjL (uu - ^ q 2 ) - = 


2V(vv - uu) - ruu ^ uv = 0 
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~~ 3V q ~~ 

- ruw t- - f vw = 0 


3r 


A 


- 2ruw |H . a (WK . i q 2) - 2bg_ . 0 


( 12 ) 


where V = V(r)/r , W = W(r)/r , uw = uw/r and vw = vw/r 
Assume the velocity correlations are of the form 

= a 2 L 3V^ 2 


u l U j 


3r 


U i U J 


(13) 


where U.U. Is a dimensionless number (turbulent correlation 
coefficient ) .Then, substituting into Eq. (12) yields a nonlinear 
system of algebraic equations for U-^Uj • The solution of this 
system is dependent upon two parameters 


N 



the swirl parameter 


3W 



the axial parameter 

(14) 


Solutions for QQ = U A Uj_ are shown in Fig. 1. Values of QQ < 0 
denote regions where turbulence is damped. With no swirl , N - 0 
and it may be shown that q 2 = 2A 2 (3W/3r) 2 and, therefore, axial 
shear is always destabilizing. This analysis 7 has been used pre- 
viously to explain the laminar central region in aircraft vortices. 
It has been shown that the swirl can in fact suppress the produc- 
tion of turbulence in the viscous core, thereby reducing the dis- 
persal of smoke which is seeded into the vortex for visualization 
purposes . 


As can be seen, the turbulent velocity correlation coefficient 
which is related to the eddy viscosity has a strong dependence on 
the parameters A and N which in turn are functions of mean flow 
gradients. It would appear unlikely that the functional form of 
an eddy viscosity could be chosen to include flow stability without 
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Figure 1 
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Swirl parameter, N 


Isopleths of constant QQ showing the dependence on 
N = ?/(r|i) and A - |£/(r §£) (Ref. 7) 


8 



resorting to a second-order closure model initially. It is to be 
expected, based on the results shown in Pig. 1, that when swirl 
is added to a jet significant changes in jet structure will result. 

Second-order turbulent modeling provides a method to compute 
the evolution of turbulent flowfields. Flow variables which are 
obtained are the mean velocity components as well as the second- 
order velocity correlations and the turbulent macroscale parameter 
A . In the next section, we will show a comparison between pre- 
dictions made by the second-order closure model and measured tur- 
bulent data. 


III. COMPARISON WITH JET DATA 


Over the years several computer codes have been developed to 
solve the second-order closure model equations for various flow 
geometries. Here, we will present results from TDV (three-dimen- 
sional vortex) a code developed by Sullivan 8 to compute the decay 
of an axisymmetric vortex. The streamwise direction is z which 
is the marching direction^ and r is the radial coordinate. We 
will use the code to predict the evolution of an axisymmetric tur- 
bulent jet and make comparisons with Wygnanski and Fiedler 9 data. 

At z = 0 the axial velocity profile was taken to be a 
Gaussian and the code was run with the axial velocity held fixed 
until the second-order velocity correlations establish themselves. 
The integral scale parameter A was taken to be constant during 
this time. The mean variables were then unlocked and solutions 
were computed to 112 diameters downstream. Figures 2-7 display 
second-order closure results compared with Wygnanski and Fiedler 
data and other theoretical results. 


The mean axial velocity distribution across the jet is shown 
in Fig. 2, at z = 42 and 57 diameters downstream of the nozzle. 
The Wygnanski and Fiedler data is presented as a composite of this 
data at z = 40, 50 and 60 diameters, since all of this data falls 
on the same solid curve. The velocity scale in this figure is the 
maximum axial velocity in each axial station and the abscissa is 
the distance from the axis of symmetry normalized by the distance 
from the nozzle. The agreement with the Wygnanski and Fielder 
data is excellent except at values of r/z > 0.15 . This error is 
due to the fact that the computational procedure requires that the 
solution asymtote to a small but finite value of W as r ■+ 00 . 

In the calculation shown in Fig. 2, this value was taken to be 
W/W m ay (z = 0) = 0.005 . A calculation with W/W ma y (z = 0) 
asymptoting to 0.02 as r 00 gave an error approximately four 
times as large for r/z > 0.15 but excellent agreement for r/z < 
0.15 . 


The variation of the fields with downstream distance is 
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shown in Pig. 3. Again, the theoretical results are in good 
agreement with the data. 

Figures 4-7 display the radial distribution of the various 
turbulent intensities at various stations downstream. In addition. 
Fig. 4 displays measurements of Corrsin and Uberoi 10 at z/D = 20 . 
Wygnanski and Fiedler attribute the 25% disagreement between their 
measurements and Corrsin and Uberois to: (1) the fact that at 

z/D = 20 the flow was not yet self -preserving, and (ii) the poorer 
response of the equipment used by Corrsin. The results of the pre- 
sent study are in significantly better agreement with those of 
Wygnanski and Fiedler in the midrange of r/z , For r/z > 0.15 , 
the major part of the disagreement is most likely due to the as- 
ymptotic condition on W as discussed above. 

In general, model predictions at least for this simple flow- 
field are in general agreement with measured data. Additional com- 
parisons may be found in Ref. 5. 


IV. COUPLING TURBULENT JET STRUCTURE 
TO AN ACOUSTIC MODEL 


We have chosen to use the RIbner and equivalent Lighthill for- 
mulation which models the sound radiated from a turbulent flow. 
Other competing models exist and our choice has been made solely 
on the grounds that the formalism behind this model is generally 
understood in the acoustic community. 

The acoustic power radiated in the 4> direction per unit 
solid angle from an elemental volume of turbulence located at y 
is 


P(4>,y) 


where 

R x U>y,f , t) = p 2 [u 2 u^. 2 + 4U x U^ u x u* + other terms] 

U x is the time mean velocity component in the x direction and 
u x is ^.he fluctuating component. The primed variables denote the 
£oint y - f/2 at time t and the unprimed denotes the point 
y + X/2 at time t + t . The omitted terms are argued to be 
either zero upon time differentiation or small. At this point, 
Ribner introduced the assumptions of isotropy and homogeneity of 
turbulence within a correlation volume in order to make further 
progress. These assumptions are not necessary, however, since 


l6ir^p a - 
o o 


/ 2j 

^-ZC R x (<J>,y,t,T) d % 


(15) 
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Figure 4. Intensity of ilw ,u fluctuation across the jet. 

, Corrsin, z/D = 20 j Wygnanski & 

Fiedler composite data z/D = 50 , 60 , 75 ; 

; present study z/D = 57 , 112 
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second-order closure modeling predicts both anisotropy and inhomo- 
geneity. At this point, however, we have chosen not. to alter the 
basic Ribner analysis, so that the acoustic model will remain simple 
It is suggested, however, that future efforts should utilize the 
full power of second-order closure. Requiring that the joint pro- 
bability density is Gaussian yields that 


2.2 2 

u x u x = U 1 


+ 2Rf 


(17) 


where Rj = u x uj^ = uquj_ . Utilizing results from the isotropic 
theory of turbulence, it was suggested that 


R 1 (1,t) 


u 2 [1 - tt5 2 /A 2 + (tt£; 2 /A 2 )] exp(-Tr^ 2 /L 2 


o>~t ) 
f (18) 


where A is the integral scale length of the turbulence and <o- 
is the inverse turbulent correlation time scale. Substituting 
Rq into (17), and (17) into (15) yields 


P(4>,y) 


. 3 4 
A iij, 
__ 

tt pa;: 
o o 




(cos 2 <f> 


+ cos 



C 

(19) 


after parrying out the time differentiation and the integration 
over £ and setting t = 0 . Also, we have approximated 


U U' = W 2 (y ) + w 
xx w 3r 


y 


4 


+ 


( 20 ) 


The convective factor 

C = 


C is given by 

2.2 

(1-5- cos4>) 2 + 

a o itaf 


-5/2 


( 21 ) 


and has been added by Ribner after the fact to account for the eddy 
igotion down the jet. The equation is yet to be integrated over 
y . This is accomplished numerically in our turbulent code. What 
remains to be done is to argue the functional form of a)~ . Second 
order closure in its present form makes no direct prediction of 
but we argue that the turbulent correlation time scale should be 
proportional to L/q where L is a characteristic turbulent scale. 
We, therefore, assume 

a) f = L/q (22) 
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Now L is taken to be proportional to our turbulent macroscale 
parameter A so that L = kA . 

Our strategy now is to determine k from data and then use 
Eq. (19) to make predictions of sound power radiated for more 
complicated flowfields. We have chosen to use Lush 11 data to make 
this calibration. Shown in Pig. 8, is a comparison of directivity 
as a function of emission angle for three jet velocities. It is 
suggested that k = 0.5 gives a best fit to the data. One does 
expect, as a consequence of the incompressible turbulence model 
and assumptions implicit in -the acoustic model, that discrepancies 
will occur as jet Mach number approaches one from below. We have, 
therefore, biased our selection of k to give better agreement at 
the lower jet velocities. Equation (19) with the second-order 
closure model will now be used to make predictions of directivity 
and total sound power radiated from more complicated jet flowfields. 

V. NOISE ESTIMATES FROM COLD JETS 

In this section the noise radiated from several cold jets are 
estimated using the models described in Sections II and III. 

Recently Lu, Ramsay and Miller 12 have made measurements of 
the noise radiated from a swirled and nonswirled model jet. We 
have attempted to numerically simulate their experiment and predict 
the total sound power radiated. In Pig. 9 is shown the measured 
swirl angle as a function of radius as well as the measured velo- 
city profile at 2 jet diameters downstream. The dashed curves are 
the predicted distributions. Unfortunately, Lu et al. have not 
presented any measurements made further downstream. To make this 
comparison we have run TDV in a pipe, so that the initial condi- 
tions at z/D* = 0 are those appropriate to a fully developed 
pipe flow. Tnen this mean velocity axial distribution was altered 
by 


W = cosa W . (r) 

pipe 

V = sina W plpe (r) (23) 

where a is the swirl angle measured at z = 0 . The predicted 
profiles are shown in Pig. 10. We also present the results of the 
calculation where the mean axial flow was not swirled in Fig. 11. 
One obvious difference is the rate of decay of the maximum 
axial velocity as a function of downstream distance. This varia- 
tion is shown plotted in Fig. 12. In general, it is also shown 
that the turbulent kinetic energy at equal downstream stations is 
higher in the jet which is swirled. This, of course, explains the 
more rapid aging of the swirled vs. the nonswirled case. We have 
estimated the thrust loss in the swirle-' case to be about 3 % . 
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Emission angle, <f> (deg) 

Figure 8. Directivity of sound intensity (dB, re 10 -12 W/m 2 ) 
(k = 0.5) vs. Lush data for three values of jet 
velocity. Lush data: single points. Simulation: 
solid line. (<J> is measured from the downstream 
jet axis.) 
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Figure 12. Predicted decay of centerline axial velocity 
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In Pig. 13, we have shown the computed directivity for the 
swirled and nonswirled jet for pressure ratio equal to 1.2 (thrust 
estimated to be 33 pounds) and 1.8 (thrust estimated to be 100 
pounds). In both cases, it is shown that the swirled jet is approx- 
imately one dB (sound intensity) noisier. At a given pressure 
ratio the thrust between swirled and nonswirled is approximately 
the same. Lu et al. have not given directivity plots for this case 
but have presented spectra 70° from the downstream axis. Here we 
estimate from their data a 5 db increase in intensity when the 
jet is swirled. To explain this result we have examined the sen- 
sitivity of our directivity prediction to small changes in turbu- 
lent initial conditions. Shown in Pig. 14 , is the computed 
directivity which results when the initial specification of the 
second-order velocity correlations are increased or decreased by 
20$. It is not difficult to argue that the swirl vanes do in fact 
increase the initial turbulence level, although an estimation of 
the amount would be difficult. Also, it seems reasonable that 
swirl vanes will alter the initial integral scale. However, this 
effect was not investigated. Within the limitations which are 
detailed above, these results do indicate that the Lu et al. cold 
swirling jet does not offer any sound intensity reduction over 
the nonrotating counterparts. 

We have also made a simulation of the jet flowfield of a high 
by-pass ratio jet engine (by-pass ratio 3 typical of a JT 15 D) 
where the primary exhaust is swirled. We have simply added a by- 
pass jet flow around the simulation of the Lu et al. model jet above. 
Here we have estimated initial conditions for the by-pass flow by 
computing the fully developed pipe flow in an annular pipe to approx- 
imate the exit conditions on the by-pass flow. In Pigs. 15 and 16 , 
is shown the initial, as well as downstream, predicted profile for 
a swirled and nonswirled primary flow. Note the initially high 
turbulence level in the region of mixing between the primary and 
by-pass flow and also the reduction of A in these regions. Again, 
we note that the effect of swirl is to increase the decay rate of 
the centerline axial velocity. Predicted directivity is shown in 
Pig. 17. Here thrust was held constant for both swirled and non- 
swirled cases at 113 pounds. Again, it is seen that swirl does not 
offer any reduction in sound power radiated over the nonswirled 
case. 


As a final computation which demonstrates the general jet flow- 
fields which can be studied using second-order closure modeling, 
we have simulated a multitube suppressor. We have used nine 
Gaussian jets with the initial turbulence specified using super- 
equilibrium theory. These initial conditions were inputted into 
a code named "WAKE" which is a nonaxisymmetric version of TDV. 

WAKE is a fully elliptic computation in the cross plane and is 
parabolic in the streamwise direction. The computed decay of axial 
velocity is shown in isopleth form in Pig. 18. Isopleths of turbu- 
lent intensity are also shown in Pig. 19. We have computed the 
directivity of this jet and compared it with that predicted from 
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Figure 13. Intensity (dB re 10 W/m ) measured from the downstream axis for 

the swirled and nonswirled jets for pressure ratios of 1.2 and 1.8 
(from Lu, et al data simulation) 
























w'/W 



Figure 15e. Predicted w'w’ profiles (no swirl) for a high 
by pass ratio jet as a function of downstream 
distance (Lu, et al simulation). See Fig. 15a 
for legend 
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Figure 17. Predicted directivity for high by pass swirled 
and nonswirled Jets, intensity (dB re 1 0~ 12 
W/m^) 
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Figure l8a. The curves are numbered such that .the maximum 

value is scaled to be 10 with 5 denoting 50 %, etc. 
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x/D = 2 .3 
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□ .s 1 1.5 2 2.5 

2 y/D 

Figure l8e. Isopleths of axial velocity at downstream position 
x/D = 9 
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an equal thrust, equal diameter axisymmetric Gaussian jet, using 
TDV. The reduction of intensity Adb is shown in Pig. 20. 


VI. THE EFFECT OF HEAT AND SWIRL ON 
TURBULENT JET STRUCTURE 


It is apparent that model predictions show that swirl does 
substantially alter the structure of a cold jet but that substan- 
tial reduction in the sound power radiated has not been computed 
to occur. It is clear that this result is opposite of that ob- 
tained by Schwartz 13 experimentally where swirl was imparted to 
hot jets. Here, he has shown that by swirling the exhaust from 
turbo-fan and turbo-jet engines, substantial reductions in total 
sound power radiated can be achieved with little loss of thrust. 

In this section we propose a possible explanation of this dis- 
crepancy based on the neglect of mean density gradients. Our ap- 
proach is to again examine the superequilibrium limit of our 
turbulent model including the effects of radial variation of den- 
sity. 


Consider the flow geometry shown in Fig. 21. To reduce the 
algebraic manipulations, we will neglect streamwise dependence on 
all variables and seek only a radial balance between terms. 
Furthermore, we will suppress the turbulence produced by the axial 
velocity and only examine the production of turbulence by swirl 
and radial density gradients. This assumption again reduces the 
algebraic complexity but still retains the swirl and radial den- 
sity gradient production terms in the rate equations for the cor- 
relations. When the Mach number is moderate and the superequili- 
brium limit of the modeled rate equations are taken, the rate 
equations become 
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Adb= 10 log P/P G 


Figure 20. Comparison of the directivity of the multitube 

suppressor to the axisymmetric jet. I = intensity 
of the multitube, Iq = intensity of the.axi- 
symmetric Gaussian 




Figure 21. Radial density, axial and swirl velocity distributions 
in a jet exhaust. 
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A,b,s are model constants which have the values 0.75 , 

0.125 and 1.8 , respectively. Lambda (A) is the integral scale 
parameter in the turbulent model and has the dimensions of length. 
Solution to the above system can be obtained by assuming 
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and substituting Eq. (25) into Eq. (24), yielding 
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The solutions to Eq. (24) are lengthy and only the 
Q 2 , UV and RU are displayed below. 
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Since Q must be positive definite, there is a relationship be- 
tween Richardson number and N for which Q 2 = 0 . On Pig. 22 
are shown Q 2 isopleths. As can be seen, self-sustained turbu- 
lence is possible only for those values of S and N which fall 
below the Q 2 = 0 curve. Above this curve, turbulence is damped 
by centrifugal effects of the swirling flow. 


With the radial density set equal to zero (S = 0) , self- 
sustained turbulence is possible if -0.68 < N < 0.18 as has been 
shown by Donaldson and Bilanin . 14 Simple parallel shear flow 
corresponds to N = 0 . Here, the critical Richardson number above 
which self-sustained turbulence is not possible is approximately 
0.55 , which agrees with that found by Lewellen, Ref. 5 . 

In Figs. 23 and 2 4 are shown curves of constant UV and RU , 
respectively. As can be seen, RU < 0 as it must be to insure 
that RR is positive definite (Eq. (26)). 


The values of S and N which may be attained in a hot 
swirling jet cover N < 0 and all S if both positive and nega- 
tive density gradients are allowed. Negative gradients might occur 
in an advanced engine which exhausts hot gas through an annular 
region. Therefore, it appears that the presence of swirl, coupled 
with mean density gradients, can profoundly affect the ability of 
the jet to produce turbulence. We, therefore, suggest that the 
mechanisms which result in noise suppression by swirling the hot 
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Figure 22. Isopleths for Q 2 (superequilibrium limit) as a 

function of swirl number N and Richardson number S . 


I 




Figure 23. Isopleths for UV (superequilibrium limit) as a 

function of swirl number N and Richardson number S 
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v/r 


(v r — v/r) 


Figure 24. Isopleths for RU (superequilibrium limit) as a 

function of swirl number N and Richardson number S 
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exhaust flow of a jet engine may be intimately related to the 
modification of the turbulent stability of the jet. 


VII. CONCLUSIONS 

A first look has been made at using second-order velocity 
correlations predicted from second-order closure turbulent trans- 
port theory to predict the generation of aerodynamic sound. It 
is suggested that second-order turbulent transport theory, in its 
present form, does contain turbulent statistical quantities which 
appear in several formulations of the aerodynamic sound generation 
problem. It is also suggested that second-order closure transport 
theory might be extended or modified to directly compute fluctua- 
tions which appear explicitly as source terms in current aero- 
dynamic sound formulations. 

Several calculations of complex jet flowfields have been 
undertaken to illustrate the versatility of second-order closure 
turbulent modeling. These computations include swirled and non- 
swirled jets and a flowfield simulating exhaust from a multitube 
nozzle. It has been shown, with second-order closure turbulent 
transport theory coupled with Ribner's formulation of the aero- 
dynamic sound generation problem, that swirling the exhaust of a 
cold jet does not offer any appreciable sound power reduction over 
a nonswirled equal thrust jet. It is suggested that the stabiliz- 
ing/destabilizing effect present when the jet flow is hot and is 
swirled, which is not included in our current model, may in fact 
explain this apparent contradiction with measured data. 
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15. Supplementary Notes 


16. Abstract 

Current formulations of the generation of aerodynamic sound by turbulence 
all require statistical information with regard to the turbulent flowfield. 
Second-order closure turbulent modeling in its present form is shown to 
provide some of this statistical information. 

The Ribnerl formulation of the generation of aerodynamic sound is coupled 
with predictions of second-order velocity correlations and integral scale 
to estimate the sound radiated from several complicated jet flows. In 
particular, it is shown that the sound radiated from a cold swirling jet 
is greater than from its nonswirling equal thrust counterpart. 

The noise radiated from the flowfield of a multitube suppressor is estimated 
and compared with an equal thrust equal diameter Gaussian jet. It is shown 
that the multitube concept is indeed quieter. 
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